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Active  and  Passive  Remote  Sensing  of  Ice 


Principal  Investigator:  Jin  Au  Kong 


Semi-Annual  Progress  Report 


This  is  a  report  on  the  progress  that  has  been  made  in  the  study  of  active  and 
passive  remote  sensing  of  ice  4*ader~thtrT»ponGorship  of  QNR  Contract  "NOOQ14-83  jfe> 
during  the  period  of  February  1,  1984  -  July  31,  1984.  During  this  period  we 
have:  (1)  calculated  the  emissivity  of  a  two-laver  anisotropic  random  medium;  and 

(2)  participated  in  the  microwave  sea  ice  measurement  program  at  the  Cold  Regions 
Research  and  Engineering  Laboratory  (CRREL). 

The  emissivity  of  a  two-layer  anisotropic  random  medium  has  been  calculated 
using  the  dyadic  Green’s  function  for  a  two-laver  anisotropic  medium  and  the  first-order 
Born  approximation.  The  emissivity  is  calculated  by  obtaining  coherent  and  incoherent 
reflectivities  and  by  making  use  of  the  relationship  e  -  l  -  r.  The  incoherent  reflectivity 
is  obtained  by  integrating  over  the  upper  hemisphere  the  bistatic  scattering  coefficients 
obtained  under  the  Born  approximation.  The  theoretical  results  are  illustrated  by 
plotting  the  emissivities  as  functions  of  viewing  angles  and  polarizations.  They  are  used 
to  interpret  the  passive  microwave  remote  sensing  data  from  vegetation  canopy  which 
also  show  strong  anisotropic  dependencies.  Data  obtained  from  the  field  measurements 
with  corn  stalks  arranged  in  various  configurations  with  perferred  azimuthal  directions 
are  successfully  interpreted  with  this  model.  It  is  expected  that  the  radiometric  data 
from  CREEL  measurements  will  also  be  successfully  interpreted  with  this  theoretical 
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model.  A  manuscript  has  been  prepared  for  submission  to  a  journal  for  publication  a 
Appendixj. 

We  have  participated  in  the  winter  microwave  remote  sensing  measurements  at 
CRREL.  Several  trips  were  made  to  the  experimental  site  at  CRREL  in  order  to  assist 
and  participate  in  the  experimental  efforts.  Our  involvement  in  the  experimental  efforts 
has  provided  us  with  valuable  insights  in  the  development  of  theoretical  models  and 
data  interpretation. 
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PASSIVE  MICROWAVE  REMOTE  SENSING  OF  AN 
ANISOTROPIC  RANDOM  MEDIUM  LAYER 
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Department  of  Electrical  Engineering  and  Computer  Science  and 
Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


Abstract 

The  principle  of  reciprocity  is  invoked  to  calculate  the  brightness  temperatures  for  passive 
microwave  remote  sensing  of  a  two-layer  anisotropic  random  medium.  The  bistatic  scatttering 
coefficients  are  first  computed  with  the  Born  approximation  and  then  integrated  over  the  upper 
hemisphere  to  be  subtracted  from  unity,  in  order  to  obtain  the  emissivity  for  the  random  medium 
layer.  The  theoretical  results  are  illustrated  by  plotting  the  emissivities  as  functions  of  viewing 
angles  and  polarizations.  They  are  used  to  interpret  remote  sensing  data  obtained  from  vegetation 
canopy  where  the  aniostropic  random  medium  model  applies.  Field  measurements  with  corn  stalks 
arranged  in  various  configurations  with  preferred  azimuthal  directions  are  successfully  interpreted 
with  this  model. 
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1.  Introduction 


In  passive  microwave  remote  sensing  of  earth  terrain,  the  random  medium  model  has  been 
developed  to  account  for  the  volume  scattering  effects  that  are  found  to  be  extremely  important. 
Both  the  wave  theory  [  1  — 31  and  the  radiative  transfer  theory  [4-5J  have  been  applied  in  the  calcula¬ 
tion  of  the  scatttering  coefficients  for  active  remote  sensing  and  the  emissivitties  for  passive  remote 
sensing. 

For  terrain  media  such  as  sea  ice  [6-7],  vegetation  fields  with  row  structures  [81,  and  vegetation 
canopy  with  preferred  azimuthal  directions  [9],  the  random  medium  must  be  characterized  with 
an  anisotropic  permittivity  tensor.  Lee  and  Kong  [10]  have  considered  anisotropy  in  both  the 
background  permittivities  and  the  randomly  fluctuating  parts.  The  backscattering  coefficients  for 
active  remote  sensing  were  calculated  and  applied  to  the  interpretation  of  radar  measurements  from 
Arctic  sea  ice. 

In  this  paper,  we  will  apply  the  principle  of  reciprocity  to  calculate  the  emissivities  [Hi  for 
passive  microwave  remote  sensing.  Making  use  of  the  dyadic  Green’s  function  for  a  two-layer 
anisotropic  medium  12j  and  the  first  order  Born  approximation,  we  first  obtain  analytical  expres¬ 
sions  for  the  bistatic  scattering  coefficients.  We  then  integrate  over  the  upper  hemisphere  and 
subtract  from  unity  to  obtain  the  emissivities. 

The  angular  responses  of  the  emissivities  are  illustrated  numerically  by  studying  their  depen¬ 
dence  on  the  variance,  the  correlation  lengths,  and  the  mean  permittivities,  all  of  which  can  be 
anisotropic.  The  theoretical  results  are  then  used  to  interpret  experimental  data  collected  from 
vegetation  fields  where  cut  corn  stalks  are  laid  on  the  ground  to  emphasize  their  azimuthal  de¬ 
pendence.  The  measured  radiometric  brightness  temperatures  will  be  shown  to  match  with  the 
theoretical  results  along  different  azimuthal  angles  of  observation 


2.  Formulation  of  the  Problem 


Consider  an  electromagnetic  plane  wave  with  a  linearly  polarized  time-harmonic  field 


E0t(f,t)  =  E(he'^'r-^ 


(2.1) 


incident  upon  a  layer  of  an  anisotropic  random  medium  with  a  permittivity  tensor 

f,(r)  =<  f,(r)  >  +?,/(r)  (2.2) 

where  <  <i(r)  >  is  the  mean  part  and  <i/(r)  represents  the  randomly  fluctuating  part  [Fig.  1). 
For  the  statistically  homogeneous  medium,  <  t  i(r)  >  will  be  a  constant  independent  of  position, 
and  1  i/(r)  is  a  random  function  of  position  whose  ensemble  average,  <  e  i/(r)  >  =  0.  It  is  assumed 
that  the  amplitude  of  (>/(r)  is  small  compared  with  <  fi(r)  >  .  In  general,  both  <  ti(r)  > 
and  tjf(r)  are  taken  to  be  uniaxial  with  the  optic  axis  tilted  off  the  z  -axis  by  angles  v  and  t pf, 
respectively,  as  shown  in  Fig.  2.  The  expressions  for  <  e  j (r)  >  and  <  <i/(r)  >  are  given  in  Ref. 
10|. 

The  layer  of  random  medium  has  boundaries  at  z  =  0  and  z  =  -d .  The  upper  region  is  free 
space  with  permittivity  e(J .  The  lower  region  is  homogeneous  and  isotropic  with  permittivity  e2- 
All  three  regions  are  assumed  to  have  the  same  permeability  n  . 

Expressing  the  electric  fields  in  all  regions  in  terms  of  the  dyadic  Green  s  functions  and  applying 
the  first-order  Born  approximation  as  explained  in  detail  in  10i,  the  first-order  scattered  electric 
field  intensity  in  region  0  is  given  by 


E'n(f)  *  >*<  E\!\f)'3  >=  jj  d*f ,d8 


rt 


{  [C0,(r,r,)  •  Q(r.)  ■  E\°\r ,)]  •  [c01(r.  f2)  •  Q(f3)  ■  E\°](r2)])  (2.3) 
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f*  I* 


where  C0i(r,  rj)  is  the  dyadic  Green’s  function  with  the  observation  point  in  region  0  and  the 
source  point  in  region  1,  Q{r)  =  u/2pe  \j(f)  is  the  fluctuation  tensor  which  stands  for  an  effective 
source  distribution,  and  E^\f)  is  the  unperturbed  electric  field  in  region  1  in  the  absence  of 
fluctuation.  C?oi(r,ri)  and  e\°\?)  are  obtained  in  Section  3  of  ,  10]  and  the  correlation  functions 
of  the  tensor  Q  are  explained  in  Section  4  of  [10] . 


3.  Bistatic  Scattering  Coefficients  and  Emissivity 


After  the  p,  p2  and  0  x  -  integrations  are  performed  as  shown  in  [lOj,  (2.3)  reduces  to 


<  Eo(r)\  >  = 


A  So,!  * 


14  r  oo  /■  r  0 

1  ?/*•//>■*■ 


(3.1) 


where  <I>  is  the  spectral  density  of  the  correlation  function  and  all  the  constants  and  functions  are 
defined  in  10  Now  we  have  to  evaluate  the  integral  of  the  following  form: 


I^f  d3z  j  dzx  [  dz2e'{-k'‘+k't*>-0*)z'e-t{-k':'k''*'-l3‘'*"-<t>(P,)e-,0A*‘-*-)  (3.2) 
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The  /?*  -integration  will  be  performed  term  by  term  using  the  contour  integration  of  residue 
calculus.  We  pick  up  the  dominant  pole  contributions  under  some  approximations  reasonable  in  the 
microwave  remote  sensing  of  earth  terrain,  as  described  in  the  paper  jlO  .  The  value  of  integration 
is  significant  only  when  p  =  s  and  q  =  t  ,  or  when  p  =  and  q  -  s'"  for  the  backscattered 
direction  The  final  result  is  the  following. 
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2  x*(/3,  =  -k» 


2n<t>{/3 


-  kn 

Z  ~  *1*1 


k'L) 

k'L) 


.2<  -*l'»  **  "if,  I'1  _ 
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for  p  =  s  .  q  —  I 

ior  p  =  ,  q  —  s^and  backscattering 

only  when  k\z  —  kjzt  and  k[z  =  k\zi 


(3.5a) 

(3.56) 

(3.5c) 


where  the  operator  +  transforms  od  ou  and  ed  e u ,  and  where  /  denotes  the  real  part  and 
n  denotes  the  imaginary  part  of  that  quantity. 

Making  use  of  (3.5)  for  evaluation  of  the  integral  in  (3.3).  we  obtain  the  bistatic  scattering 
coefficients,  which  are  defined  by  Peake  11  as 


_  _  4ir r2  <  E0(r)  >a 

lja(kf,k0 ,)  =  lim  - — — 2 — 

A  cos  9 o%  Eqi  3 


•j  .  ....  2 

where  L',„  4  is  the  incident  electric  field  intensity  with  polarization  J.  <  E0(r)  is  the 

mean  scattered  field  intensity  with  polarization  a.  k,„  is  the  incident  wave  vector  at  incidence 
angles  (0.„,o,„)  and  k,  is  the  scattered  wave  vector  at  scat tered  angles  (0f.of)  .  The  final  results 
are  as  follows: 


o3.l« 
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(3.6) 
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i  1 


(3.7) 
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where  the  expressions  for  J,  s  are  given  in  the  Appendix. 

By  using  energy  conservation  and  reciprocity  arguments,  we  can  express  the  emissivity  at  angle 
(0ot,<Pot)  with  0  -polarization  in  the  following  form  [11]: 


£  *  $01)  —  1  A p (0 ot ,  <p 0i)  .  0  —  H  or  V 


(3.8a) 


—  1  r9c  (@ot  •  fiat )  r  di  (^<>«  -  *Po\ ) 


(3.86) 


where  A0(9ot,  <pot)  is  an  “albedo"  defined  as  the  fraction  of  the  power  incident  on  the  surface 
from  the  direction  (0,„,  <£,>,)  at  3  -polarization  that  is  reflected  or  rescattered.  r/?c(0in,  <p„,)  is  the 
"coherent  reflectivity”  which  denotes  the  fraction  of  the  power  reflected  in  a  specular  direction,  and 
is  the  “incoherent  reflectivity”  which  denotes  the  fraction  of  the  power  scattered  over 
the  upper  hemisphere.  The  reflectivities  are  given  by 


r0c{9oi  ,  <Pol  ) 
r3x(0Oi,<i>m) 


^  (^ot  j  <?ot )! 

a  =  H,V 
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|  /*2*  ri 
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where  ~j0a  is  given  in  (3.7)  and  B0at  is  calculated  in  Ref  12  as  follows: 


(3.9a) 

(3.96) 


B  i-iry  ,  B\  1 1  /l  o  “**  -A  0  o  (  ^  1  O  0  2^  (>)  ^  '(Of  Xao(L,Of  Xf’2^  r  )  f  u 

~X0e{L2Ol,  -  XUe0)X0a  -  X0t(L>o,  -  M\tt)Xta  (3.10) 


where  all  the  coefficients  are  defined  in  12  and  should  be  evaluated  at  angles  of  incidence  (0O,.  <£,„) . 


4.  Discussion  and  Application  of  the  Results 

First,  it  can  be  easily  shown  that  (3.7)  reduces  to  (5.24)-(5.32)  of  the  paper  }  10!  when  9  f  =  9ol 
and  O.  -  7r  -  <6„,  .  which  are  the  results  for  the  backscattering  coefficients.  Also,  in  the  limiting 
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case  when  both  the  mean  and  fluctuation  are  isotropic,  i.e.  when  ti*  =  ti,  t\zf  —  fj/,  and 
rp  =  tpf  —  0 ,  the  result  (3.7)  reduces  to 
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where 


3,  =  (kpl  -  *,)  z(k\zt  -  Ar',J 

=  (kpi  -  £,)  -  i(*i„  -  At;  j 

30  =  2k  pt 

and  A  is  defined  in  the  appendix. 

These  results  are  exactly  the  same  as  those  for  bistatic  scattering  coefficients  obtained  in  the 
case  of  two-layer  isotropic  random  media  by  Zuniga  and  Kong  [13;-  It  is  interesting  to  note  that 
our  results  can  be  shown,  with  some  algebraic  manipulation,  to  satisfy  the  following  four  relations: 


cos^„,7^Q(t,  s)  =  cos  Of-iag(s.  1) ,  3,  a  —  H  or  V 


(4.2) 


where  t. fie  first  lettt.er  in  the  brackets  indicates  the  direction  from  which  the  radiation  is  incident,  and 
the  second  letter  indicates  the  direction  into  which  the  radiation  is  scattered.  This  is  a  consequence 
of  the  reciprocity  theorem  1 1  as  the  anisotropic  random  medium  is  reciprocal  and  its  permittivity 
tensor  is  symmetric. 

As  an  illustration  of  our  theoretical  results,  the  emissivities  e/y  and  ty  are  plotted  as  a 
function  of  the  incidence  angle  at  10GHz  in  Fig.  3-4.  The  correlation  function  is  assumed  to  be 
exponential  in  both  the  lateral  and  the  vertical  directions  and  it  is  assumed  that  U-  —  n>/  as  in 
ID  The  dependence  of  the  emissivity  on  the  azimuthal  incidence  angle  <pot  as  shown  in  Fig.  3  is 
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a  clear  indication  of  the  medium  anisotropy.  When  61  >  6%  ,  the  strength  of  fluctuation  is  stronger 
in  the  lateral  direction  (the  plane  perpendicular  to  the  optic  axis)  than  in  the  vertical  direction 
(along  the  optic  axis  of  fluctuation).  Then  for  large  azimuthal  incidence  angles  ( <p0l  —  80c), 
the  H  -polarized  wave,  whose  electric  field  is  almost  directed  laterally,  is  scattered  more,  resulting 
in  smaller  emissivity.  On  the  other  hand,  the  V  -polarized  wave  is  scattered  less  for  large  0ot  , 
resulting  in  higher  emissivity. 

In  Fig.  4,  for  large  80s).  when  c'/2  is  increased  from  0.02  to  0.1,  e H  remains  ap¬ 

proximately  the  same,  but  ev  increases  noticably.  It  is  due  to  the  anisotropy  of  mean  permittivity. 
When  the  loss  in  the  vertical  direction  is  increased,  the  H  -polarization  is  not  significantly  affected, 
but  the  V  -polarization,  which  has  components  in  the  vertical  direction,  will  experience  more  of 
the  absorption  effect,  resulting  in  higher  emission. 

As  an  application  of  our  theoretical  results,  we  choose  to  interpret  the  field  data  of  series  3 
obtained  by  NASA  and  USDA  at  BARC  (Betsville  Agricultural  Research  Center)  with  their  1982 
corn  destruction  experiments  ,14-15  The  corn  stalks  above  the  ground  were  cut  and  removed,  and 
the  brightness  temperature  is  measured  with  a  radiometer  for  the  stubble  ground.  Then  the  stalks 
were  laid  down  on  the  ground  and  aligned  in  specific  orientations.  The  measurements  were  made 
with  the  radiometer  looking  parallel  and  perpendicular  to  the  row  direction  of  corn  stalks  at  two 
frequency  bands,  the  L-band  (1,4 GHz)  and  the  C-band  (5 GHz) .  Because  the  corn  stalks  have  a 
preferred  orientation  in  the  horizontal  plane,  they  are  modelled  as  an  anisotropic  random  layer  with 
the  optic  axis  in  the  horizontal  plane  (il%  =  v/  -  90c)  .  As  is  clearly  seen  in  Figs.  5-10,  the  theory 
explains  fairly  well  the  data  behavior  at  different  angles  of  incidence,  polarizations,  frequencies, 
and  orientations.  The  letters  V  and  H  denote  the  experimental  data  and  the  curves  denote  the 
theoretical  results.  All  six  sets  of  data  are  matched  with  a  single  set  of  parameters:  6,  =  0.03, 
62  =  0.075 ,  /p  =  2cm  ,  lz  —  15cm  ,  and  with  the  permittivities  of  the  stalk  layer  and  the  ground 
following  the  frequency-dependence,  as  given  by  Ulaby  and  Jedlicka  16*. 
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In  the  calculation  of  the  brightness  temperature,  the  following  relation  is  used,  assuming 
uniform  temperature  distribution  T0  ■ 


T[b0)  =  e0To  (4.3) 

where  (3  is  the  state  of  polarization  to  be  considered  and  T„  is  the  physical  temperature  of  the 
ground  in  ZK  . 

We  first  observe,  by  comparing  Fig.  6  and  Fig.  9  (stalks-covered)  with  Fig.  5  and  Fig.  8 
(stubble),  respectively,  that  the  brightness  temperatures  have  greatly  increased,  on  account  of  the 
com  stalk  layer.  Second,  by  comparing  Fig.  7  [<j>ot  =  0°)  with  Fig.  6  (<j>ol  =  90°)  at  1.4GHz, 
we  see  that  the  horizontal  polarization  increases  and  the  vertical  polarization  decreases,  both  so 
substantially  that  there  is  a  crossover  between  curves  of  TgH ^  and  .  This  phenomenon  can 

be  explained  by  the  anisotropy  of  the  dielectric  loss  of  the  com  layer  where  ("z  >  r".  For  an 
azimuthal  observation  angle  perpendicular  to  the  stalk  orientation  (<&OI  =  0C) .  the  horizontal 
polarization  has  an  E  field  parallel  to  the  stalk  and  the  vertical  polarization  has  E  perpencicular 
to  the  stalk.  Then  the  H  -polarization  experiences  more  absorptions  than  the  V  -polarization, 
resulting  in  an  increase  of  the  brightness  termperature  for  the  H  -polarization  and  a  decrease  for 
the  V  -polarization. 

Third,  we  compare  Fig.  10  ( 0OX  ~  90c)  with  Fig.  9  (c>,„  =  O'-)  at  5GHz.  It  is  seen  that 
the  horizontal  polarization  increases,  wheras  the  vertical  polarization  decreases.  The  responses  are 
completely  opposite  from  those  in  the  case  of  the  L-band  (/  =  1.4  GHz).  For  the  C-band,  where 
the  dielectric  loss  is  chosen  to  be  smaller,  the  scattering  effect  is  dominant  over  the  absorption 
effect  for  both  polarizations.  With  this  in  mind,  the  opposite  responses  can  be  explained  by  the 
anisotropic  fluctuation  of  the  stalk  layer  where  <  62  ■  For  observation  angles  parallel  to  the  stalk 
orientation  ( <pox  =  902) .  the  horizontal  polarization  has  an  E  field  perpendicular  to  the  stalk. 


It) 


whereas  the  vertical  polarization  has  E  -components  parallel  to  the  stalk.  Due  to  the  stronger 
fluctuation  along  the  stalk  (62) ,  the  H  -polarization  experiences  less  scattering  whereas  the  V  - 
polarization  experiences  more;  hence,  the  result  is  an  increase  in  Tb  for  the  H  -polarization  and 
a  decrease  for  the  V  -polarization. 


5.  Conclusions 

With  the  use  of  the  dyadic  Green’s  functions,  a  two-layer  anisotropic  random  medium  has 
been  studied  and  applied  to  the  passive  remote  sensing  of  earth  terrain  media.  Applying  the  first- 
order  Born  approximation,  the  bistatic  scattering  coefficients  are  obtained  in  analytical  forms.  The 
emissivities  are  then  computed  by  invoking  the  principle  of  reciprocity.  It  is  shown  that  both  the 
absorptive  and  randomly  fluctuating  properties  of  the  anisotropic  medium  affect  the  behavior  of 
the  resulting  brightness  temperatures  both  in  theory  and  in  actual  controlled  field  measurements. 
The  theoretical  results  are  favorably  matched  with  the  experimental  data  obtained  from  the  corn 
stalks  with  detailed  ground  truth  information.  When  and  if  such  ground  truth  is  available  for  sea 
ice  measurements,  the  theoretical  model  is  ideal  for  the  interpretation  of  the  experimental  data  as 
sea  ice  is  also  a  highly  anisotropic  random  medium. 
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Scattering  geometry  of  a  two-layer  anisotropic  random  medium. 

Geometrical  configuration  characterizing  a  permittivity  tensor. 

Incidence  angle  responses  of  emissivity  due  to  the  effect  of  anisotropic  permit¬ 
tivity  fluctuation  when  >  S2  ■ 

Incidence  angle  responses  of  emissivity  due  to  the  effect  of  anisotropic  back¬ 
ground  permittivity  when  <j>ot  =  80°  . 

Interpretation  of  radiometric  data  for  ground  with  stubbles  at  L  band  . 
Interpretation  of  radiometric  data  for  cut  stalks  parallel  at  L  band  . 
Interpretation  of  radiometric  data  for  cut  stalks  perpendicular  at  L  band  . 
Interpretation  of  radiometric  data  for  ground  with  stubbles  at  C  band  . 
Interpretation  of  radiometric  data  for  cut  stalks  perpendicular  at  C  band  . 


Figure  10  Interpretation  of  radiometric  data  for  cut  stalks  parallel  at  C  band 
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